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Abstract. This contribution briefly reviews the reverberation mapping tech- 
nique which leads to determination of black hole masses. I focus on the emerging 
relation between the broad-line region size and the active galactic nucleus (AGN) 
luminosity, and on an overview of recent results of reverberation mapping stud- 
ies which are starting to cover the full AGN luminosity range. Preliminary 
results and time lag determination from a reverberation mapping program of 
high-luminosity quasars are also presented. 



1. Introduction 

While the physical origin of Active Galactic Nuclei (AGNs) continuum variability 
is still unclear, it is possible to use the "reverberation mapping" technique to 
study the geometry and kinematics of the gas in the Broad Line Region (BLR) , 
and to deduce the mass of the central black hole (BH) in the center of the 
AGN. The technique is based on the response of the ambient gas to changes in 
the central continuum source. Such response was first used to account for the 
observations of the apparent expansion of Nova Persei in 1901 (Couderc 1939) 
and has been proposed to explain the light curves of Type I supernovae (Morrison 
& Sartori 1969). The method was first suggested to be used in the analysis of 
AGN light curves by Bahcall, Kozlovsky, & Salpeter (1972), who calculated the 
response of the line intensity in a spherical distribution of gas. Blandford & 
McKee (1982) were the first to coin the term "reverberation mapping" and put 
it into mathematical formalism with the fundamental equation that relates the 
emission-line and continuum light curves, L(v,t) and C(t): 

L( V ,t) = f y(v,t-T)C(T)dT, (1) 



where v is the velocity field of the BLR (which manifest itself in the emission-line 
profile), and ^f(v,r) is defined from this equation as the transfer function, which 
holds in it the information about the geometry and kinematics of the BLR. The 
latter was studied by several authors who showed how ^>(v,t) can be derived 
from the observed continuum and line light curves and how ^f(v, r) will changed 
for different geometries and kinematics of the BLR (e.g., Welsh & Home 1991; 
Perez et al. 1992; Home et al. 2004). 

In practice, in order to get the line and continuum light curves, the AGN 
needs to be monitored frequently over a period of time (from days to weeks de- 
pending on the AGN luminosity and variability characteristics). First attempts 
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to carry out AGN reverberation mapping used poorly sampled light curves and 
low resolution spectra (e.g., Peterson 1988 and references therein), thus leading 
to a collapse of the two-dimension transfer function, ^f(v, r), into one-dimension 
transfer function ^(t). In fact this is further collapsed into only one parame- 
ter: the time lag between the line light curve and the continuum light curve. 
The time lag is defined as the centroid of the cross correlation function (CCF) 
between the continuum and line light curves: 

Fccf{t) = {Nucol)- 1 ]T C(t)L(t + r), (2) 

t 

where N is the number of points used in the sum for the lag r, uc and <jl are 
the rms of the light curves, and the light curves have zero mean. The centroid 
of this CCF is taken to be a measure for the size of the BLR, denotes as Rblr- 
Once Rblr is found from reverberation mapping the mass of the BH in the 
center of the AGN can be estimated using: 

M BH = JG^RblrV 2 , (3) 

where V is a measure of the BLR clouds' velocity, and / is a dimensionless factor 
that depends of the geometry and kinematics of the BLR. The mass — luminosity 
relation is further discussed in Peterson's contribution to these proceedings and 
will not be discussed here. 

Over the past two decades many monitoring campaigns were carried out 
and enabled the measurement of Rblr m about three dozen AGNs. Some of 
the notable projects are: (1) Individual monitoring of Seyfert I galaxies (e.g., 
Mrk279, NGC5548, NGC4151 — Maoz et al. 1991 — and many more by the 
'AGN Watch"0 projects — Peterson 1999). (2) The Lover of Active Galaxies 
(LAG) campaign (e.g., Robinson 1994). (3) The Ohio State University moni- 
toring program (Peterson et al. 1998). (4) The Wise Observatory and Steward 
Observatory 17 Palomar-Green (PG) quasars monitoring program by Kaspi et 
al. (2000). For recent reviews of the subtleties of the reverberation mapping 
technique see Peterson (1993), Netzer & Peterson (1997), Peterson (2006) and 
references therein. In the following sections I will summarize the current situa- 
tion and recent studies using reverberation mapping. I will also try to point on 
directions I think reverberation mapping should take in the near future. 

2. Size — Luminosity Relation 

Peterson et al. (2004) compiled all available reverberation-mapping data, ob- 
tained up to then, and analyzed them in a uniform and self-consistent way to 
improve the determination of the time lags and their uncertainties and derived 
Rblr for all objects with available data. Kaspi et al. (2005) used these size 
measurements to study the relation between Rblr and the Balmer emission line, 
X-ray, UV, and optical continuum luminosities. This relation is a fundamental 
relation in AGNs study since both quantities (Rblr and L) are directly obtained 
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Figure 1. Balmcr-line BLR size plotted versus the XL \ (5100 A) luminosity 
(in units of ergss -1 ). The BLR size of each data set is determined from the 
averaged Balmer-line time lags. Objects with multiple data sets have been 
averaged to one point per object. See Kaspi et al. (2005) for further details. 



from measurements with minimum assumptions and models. Once this relation 
is determined for the objects with reverberation-mapping data, it is used to 
estimate the mass of the BH in other AGNs by using a 'single-epoch measure- 
ments' of their luminosity and the line width (e.g., Wang & Lu 2001; Woo & 
Urry 2002; Grupe & Mathur 2004). To study the robustness of the correlation 
Kaspi et al. (2005) used data subsamples and two different regression methods: 
(1) The linear regression method of Press et al. (1992), in which a straight-line 
is fitted to the data with errors in both coordinates (known as FITEXY) and 
follow Tremaine et al. (2002) procedure to account for the intrinsic scatter in 
the relation. (2) The bivariate correlated errors and intrinsic scatter (BCES) 
regression method of Akritas & Bershady (1996). 

Assuming a power-law relation Rblr °c L a Kaspi et al. (2005) find that 
the mean best-fitting a is about 0.67 ± 0.05 for the optical continuum and the 
broad H/3 luminosity, about 0.56 ± 0.05 for the UV continuum luminosity, and 
about 0.70 ±0.14 for the X-ray luminosity. They also find an intrinsic scatter of 
~ 40% in these relations. In Fig. 1 the mean Balmer-line BLR size versus the 
\L\(5W0 A) luminosity is plotted, with one averaged point per object. Four fits 
are shown: using all 35 points, excluding the two low luminosity AGNs (which 
are heavily influenced by intrinsic reddening), and with the two fitting methods. 
Within the luminosity range of the measurements (10 43 -10 46 ergss -1 ) all fits 
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are consistent with each other and all are well within the scatter of the points 
in the plot. The disagreement of these results with the theoretical expected 
slope of 0.5 indicates that the simple assumption of all AGNs having on average 
same ionization parameter, BLR density, column density, and ionizing spectral 
energy distribution, is not valid and there is likely some evolution of a few of 
these characteristics along the luminosity scale. 

Bentz et al. (2006) used high-resolution images of the central region of 
14 of the reverberation-mapped AGNs (mostly the low luminosity ones) and 
accounted for the host-galaxy star light contamination of the AGN luminosities. 
Removing the star-light contribution and excluding several points (some of which 
do not have measured star-light contributions, do not have reliable H/3 BLR 
size measurement, or have nuclear structure and reddening that influence the 
luminosity measurement), they find the power-law slope of the size — luminosity 
relation to be 0.518 ± 0.039 (see Bentz' contribution in these proceedings). 

3. Expanding the Luminosity Range 

Current reverberation mapping studies cover the luminosity range of ~ 10 42 - 
10 46 ergs -1 . Since the full AGN luminosity range is 4 orders of magnitude larger 
than this range and span the range 

10 40_ 1Q 48 erggg-i^ there 

is an essential need 

to carry out reverberation mapping studies for lower- and higher-luminosity 
AGNs. Hopefully such broadening of the luminosity range will help to define 
better the slope of the size — luminosity relation. 

Since most reverberation-mapping studies mentioned above were based on 
Balmer emission lines, generally H/3, and on optical luminosity, single-epoch 
estimates for objects at redshifts z ^ 0.6 have had to rely either on IR ob- 
servations (e.g., Shemmer et al. 2004) or on attempts to extend the optically 
based size — luminosity relation to UV luminosities and UV broad emission lines 
(e.g., McLure & Jarvis 2002; Vestergaard &: Peterson 2006). While important 
progress has been made, there are still a number of potential problems that need 
to be addressed (e.g., Maoz 2002; Baskin & Laor 2005). Thus, the single-epoch 
measurements depend on the untested assumption that these extrapolations are 
valid. Although o posteriori explanations of the physical plausibility of the ob- 
served relations can be found, it is quite possible that subtle or strong deviations 
from the relations occur at high luminosities or redshifts (Netzer 2003). 

3.1. Reverberation Mapping of High-Luminosity AGNs 

Reverberation mapping of high-luminosity quasars is an ambitious task. Quasars 
of the highest luminosities (with bolometric luminosity, Lbol ~ 10 47 -10 48 ergs s _1 ) 
are expected to harbor some of the most massive BHs known, with Mbh ^ 
10 9 M Q . More massive BHs may have slower continuum flux variations with 
smaller amplitudes (e.g., Giveon et al. 1999; Vanden Berk et al. 2004). The re- 
quired observing periods of high-luminosity quasars are also significantly length- 
ened by cosmological time dilation, since such sources are typically found at high 
redshifts (z ^ 1). On the other hand, the ability to monitor high-z objects in the 
rest-frame UV, in which AGN variability amplitudes are routinely higher than 
in the optical, can lead to better characterized continuum light curves. The 
smaller intrinsic variability amplitude of the continuum could result in smaller 
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flux- variability amplitudes for the emission lines, affecting the ability to detect 
the time delay in the BLR response. Furthermore, high-redshift sources are 
fainter and hence more difficult to observe. Probably due to all of these possible 
problems, no reverberation measurements exist for AGNs with L ^> 10 46 ergs s _1 , 
and several attempts at such measurements have so far not been successful (e.g., 
Welsh et al. 2000; Trevese et al. 2006; A. Marconi 2005, private communica- 
tion). 

In view of the many unknowns and the opposing effects entering the above 
discussion, and considering the importance of the subject, over a decade ago, we 
began a reverberation-mapping program aimed at high-luminosity, high-redshift 
AGNs (Kaspi et al. 2007). The sample of 11 high-luminosity quasars was 
selected in 1994 from the Veron-Cetty & Veron (1993) catalog. These are high 
declination (S > 60°) objects, with observed magnitude V <; 18, redshifts in the 
range 2 < z < 3.4, and in the luminosity range of 10 46 4 <; AL A (5100 A) <; 10 47 ' 6 
ergss" 1 . This is an order of magnitude higher than other AGNs with existing 
reverberation measurements (see Fig. 1). 

All 11 AGNs are monitored photometrically at the Wise Observatory (WO) 
1 m telescope since 1995 in B and R bands. Since the targets have high decli- 
nations, they can be observed from the WO for about 10 months a year, with 
observations scheduled about once every month. Spectrophotometric monitoring 
of six of the 11 quasars has been carried out since 1999 at the 9 m Hobby-Ebcrly 
Telescope (HET). Observations are carried out using a comparison star which is 
observed in the slit simultaneously with the quasar and serves for spectrophoto- 
metric calibration of the quasar under non-photometric conditions. For details 
of the observational and reduction technique see, e.g., Maoz et al. (1994) and 
Kaspi et al. (2000). 

The continuum light curves of all 11 quasars show variations of 10-70% 
measured relative to the minimum flux. Comparing the variability character- 
istics of this sample to these of the lower-redshift PG quasar sample (Kaspi et 
al. 2000) the variability of the later is about double the variability of the high- 
luminosity quasars. The lower rest-frame variability measured in the continuum 
for the current sample is probably a manifestation of the general trend that 
high-luminosity AGNs have longer variability timescales (e.g., Vanden Berk et 
al. 2004), perhaps as a result of their higher BH masses. 

None of the four Lya light curves shows significant variability. In contrast, 
the two Ciii]A1909 light curves and all six CivA1550 light curves show signifi- 
cant variability. The variability measures of these emission-line light curves are 
comparable to, or even greater than, those of their corresponding continuum 
light curves. There are few previous AGN UV data sets with which to compare 
these possible trends. The only quasar with UV variability data of similar qual- 
ity is 3C273. Interestingly, Ulrich et al. (1993) noted the non-variability of Lya 
in this object, at a level of < 5%, over a period of 15 years, despite factor-of-two 
variations in the continuum during the same period (there are no data for the 
C IV and C III] lines during that time) . 

The main objective of our program is to detect and measure a time delay be- 
tween the continuum and the line-flux variations in high- luminosity AGNs. The 
significant continuum and line variations that was observed during a decade 
demonstrate that, at least in principle, such a measurement may be feasible. 
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Figure 2. Optical continuum and C iv light curves for S5 0836+71. Squares 
are spectrophotometric data from the HET. Triangles are photometric data 
from WO. Time is given in Julian Day (bottom) and UT date (top). Con- 
tinuum flux densities, f\, are given in units of 10~ 16 ergscm" 2 s" 1 A -1 and 
emission-line fluxes are given in units of 10" 14 ergs cm" 2 s _1 . 



Examining the light curves of the six quasars with emission-line data at the 
current stage of our project, all but one currently suffer from either low vari- 
ability amplitude in the emission-line light curves or monotonically increasing 
or decreasing continuum light curves. The one current exception is S5 0836+71; 
although the data for this quasar are still not ideal for reverberation mapping 
they do allow a preliminary measurement of the emission-line to continuum lag. 
Figure 2 shows the continuum and C IV emission line light curves for this object. 
These light curves has largest variation among all our monitored quasars. Fig- 
ure 3 shows the CCFs for these two light curves. The tentative time lag between 
the Civ line and the continuum of S5 0836+71 is found to be 595jzff days, or 
188^37 days in the quasar rest frame. 

The mean FWHM of the C IV line measured from the mean spectrum of 
S5 0836+71 is about 9700 km s" 1 . Using Eq. 5 of Kaspi et al. (2000) and the time 
lag of 188 days, the central mass of S5 0836+71 is estimated to be ~ 2.6 x 1O 9 M . 
This is the highest mass directly measured for a BH in an AGN using rever- 
beration mapping. 3C273 (=PG 1226+023), the quasar with the highest di- 
rectly measured mass so far, has a mass of 8.9 x 1O 8 M , AL A (1350 A)=2.0x 10 46 
ergss" 1 , and AL A (5100 A)=9.1xl0 45 ergss" 1 . Thus, S50836+71 has a factor 
three higher mass and a factor ~ 6 higher luminosity than 3C 273. 

3.2. Reverberation Mapping of low-Luminosity AGNs 

Reverberation mapping of low-luminosity AGNs might be considered fairly easy 
due to the short timescales involved and the expected high-amplitude variability. 
However, as these objects are of low luminosity there will be a need of a 3-10 m 
class telescope to carry out reverberation mapping campaigns for such objects 
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Figure 3. Cross-correlation functions, ICCF (solid curve, White & Peterson 
1994; Gaskell 1994) and ZDCF (circles with error bars, Alexander 1997), 
between the continuum and the C iv emission-line light curves of S5 0836+71 
from Fig. 2. 

(for example the candidate AGNs with intermediate-mass BH from the sample 
by Green & Ho 2004). So far no reverberation mapping campaign for AGNs 
with optical luminosity ^ 10 42 ergss _1 were carried out successfully. 

The one exception is NGC4395 (AL A (5100 A)=5.9xl0 39 ergss" 1 ) in which 
Peterson et al. (2005) measured the BLR size of its C iv emission line to be 
1 ± 0.3 light hr. This is consistent with the size expected from extrapolating 
the Rblr-L relation to lower luminosities. However, two optical campaigns to 
determine the H/3 time lag of NGC4395 were so far unsuccessful due to bad 
weather (Desroches et al. 2006 and contribution in these proceedings; Kaspi et 
al. 2007, in preparation). 

Until recently, only four AGNs had measured C iv reverberation time lags: 
NGC 3783, NGC 5548, NGC 7469, and 3C 390.3 (see Peterson et al. 2004, for a 
summary). NGC 4395 is four orders of magnitude lower in luminosity than those 
four AGNs and S5 0836+71 is 3 orders of magnitude higher. Thus, a preliminary 
C iv-size — UV-luminosity relation over 7 orders of magnitude in luminosity can 
be determined. Figure 4 shows the data for the above 6 objects with the best 
fit slope using the two different fitting methods; 0.551 ± 0.053 and 0.559 ± 0.025 
for the FITEXY and BCES methods, respectively. 

4. Further Prospects of Reverberation Mapping 
4.1. Two-Dimension Reverberation Mapping 

The initial goal of the reverberation-mapping technique (Blandford Sz McKee 
1982) is to study the geometry and kinematics of the BLR. This was not yet 
achieved due to poor spectroscopic data and insufficient time sampling. More de- 
tailed information on BLR geometry and kinematics in AGN can be obtained by 
studying line profile variations. Various researchers computed two-dimensional 
echo images for specific gas motions in the BLR: outward moving gas clouds, 
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Figure 4. Size-luminosity relationship based on the C IV emission line and 
the UV continuum. Linear fits to the data are shown. 

inward falling clouds, circulating gas clouds in a plane, or clouds orbiting in 
randomly-inclined orbits (e.g., Welsh & Home 1991; Perez et al. 1992; Home et 
al. 2004). 

Only few studies attempted the two-dimensional reverberation mapping 
on actual data and none resulted with conclusive results. Recently, Kollatschny 
(2003; and contribution in these proceedings) used the monitoring data of Mrk 110 
to study the variations in the line profiles and produced time-delay versus ve- 
locity maps. These maps resemble a disk transfer-function maps. Kollatschny 
(2003) finds that the outer line wing respond before the inner line profile and 
confirms the stratification of the ionizing structure in the BLR by showing that 
lower ionization lines respond after the higher ionization lines. Thus, With re- 
cent improved sensitivity of optical telescopes and the accumulated experience 
from previous mapping campaigns, it is becoming possible to obtain the crucial 
information about the geometry and kinematics of the BLR gas. 

4.2. Dust Reverberation Mapping 

IR emission in AGNs is considered to be from the alleged torus region which is 
at distances larger than the BLR from the BH. Thus, IR reverberation mapping 
might reveal the distance of the torus from the BH. Only few IR monitoring 
campaigns were carried out in the past two decades (e.g., Clavel et al. 1989; 
Glass 1992; Sitko et al. 1993). Recently, Suganuma et al. (2006; and see 
contribution in these proceedings) monitored the optical and IR emission in four 
additional objects and determined time lags. Together with previous results they 
are able to construct the torus-size — luminosity relation for 10 objects. They 
find the torus size to strongly correlate with the optical luminosity (the time lag 
is consistent with the square root of the luminosity) and that it weakly correlates 
with the mass of the BH. 
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4.3. X-ray FeKa Reverberation Mapping 

Several studies suggested the application of the reverberation mapping technique 
to the broad 6.4 keV FeKa line seen in the X-ray band (e.g., Reynolds et al. 
1999). This line is considered to emerge from the accretion disk in the very close 
vicinity of the BH, and using reverberation mapping will allow the measurement 
of the disk's size. So far several attempts to apply this method did not produce 
significant results (e.g., Ballantyne et al. 2005 and references therein), implying 
either on the complicated connection between the broad 6.4 keV FeKa line flux 
and the X-ray continuum, or on the fact that the X-ray data obtained so far 
were not sufficient for reverberation mapping. 

5. Summary 

Over the past two decades reverberation mapping of AGNs have yield mea- 
surements of the BLR size in about three dozen AGNs in the luminosity range 
~ 10 42 -10 ergs -1 ). This enables to establish a scaling relation between the 
BLR size and luminosity in AGNs which, in turn, allows the estimate of the BH 
mass in AGNs. Using reverberation mapping of different emission lines implies 
about the radial ionization stratification of the BLR (higher ionized specious 
emits from inner BLR), and that motion of the gas in the BLR are virial and 
primarily orbital. Current BLR studies should aim at broadening the luminosity 
range to all AGNs (~ 10 40 -10 40 ergs -1 ) and first steps toward low- and high- 
luminosity AGNs are being taken. Two-dimensional reverberation mapping is 
a promising direction which will produce information about the geometry and 
kinematics of the BLR. Reverberation mapping in the IR enables measurement 
of the dusty region in AGNs (torus) which seems to surrounds the BLR. On 
the other hand reverberation mapping of the inner accretion disk, using X-ray 
observations, is still to be proven feasible. 
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